Explosive Dark Matter Annihilation 
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If the Dark Matter (DM) in the Universe has interactions with the standard-model particle, the 
pair annihilation may give the imprints in the cosmic ray. In this paper we study the pair annihilation 
processes of the DM, which is neutral, however has the electroweak (EW) gauge non-singlet. In this 
estimation the non-relativistic (NR) effective theory in the EW sector is a suitable technique. We 
find that if the DM mass is larger than about 1 TeV, the attractive Yukawa potentials induced 
by the EW gauge bosons have significant effects on the DM annihilation processes, and the cross 
sections may be enhanced by several orders of magnitude, due to the zero energy resonance under 
the potentials. Especially, the annihilation to two 7's might have a comparable cross section to 
other tree-level processes, while the cross section under the conventional calculation is suppressed 
by a loop factor. We also discuss future sensitivities to the 7 ray from the galactic center by the 
GLAST satellite detector and the Air Cerenkov Telescope (ACT) arrays. 
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Nature of the Dark Matter (DM) in the Universe is an 
important problem in both particle physics and cosmol- 
ogy. The Weakly-Interacting Massive Particle (WIMP), 
X°, is a good candidate for the DM. It works as the cold 
dark matter in the structure formation in the Universe. 
High resolution iV-body simulations show that the cold 
dark matter hypothesis explains well the structure larger 
than about 1 Mpc [1]. Also, the WMAP measured the 
cosmological abundance precisely as VLqm = 0.27 ± 0.04 
[2] . Now we know the gravitational property of the DM 
in the structure formation and the abundance and dis- 
tribution in the cosmological scale. The next questions 
are the constituent of the DM and the distribution in the 
galactic scale. 

If the DM is SU(2) L non-singlet, a pair of the DM 
could annihilate into the standard-model (SM) particles 
with significant cross sections [3]. We call such DM's as 
electroweak-interacting massive particle (EWIMP) DM 
in this paper. The detection of exotic cosmic ray fluxes, 
such as positron, anti-proton and 7 ray, may be a fea- 
sible technique to search for the DM's. Since some 
DM candidates in the supersymmetric (SUSY) models 
have interactions with the SM particles, these annihila- 
tion processes are extensively studied. Especially, excess 
of monochromatic 7 ray due to the pair annihilation is 
a robust signal if observed, because the diffused 7-ray 
background must have a continuous energy spectrum [4] . 
Searches for the exotic 7 ray from the galactic center, the 
galactic halo, or even from extra galaxies are ones of the 
projects in the GLAST satellite detector and the big Air 
Cerenkov Telescope (ACT) arrays such as CANGAROO 
III, HESS, MAGIC and VERITAS. 

In the previous estimates, the cross sections for the 
EWIMP are evaluated at the leading order in the per- 
turbation. However, the DM is non-relativistic (NR) in 
the current Universe. In this case, if the EWIMP mass 
m is much heavier than EW scale, the EWIMP wave 
function may be deformed under the Yukawa potentials 



induced by the EW gauge boson exchanges and it may 
give a non-negligible effect in the annihilation processes. 
Furthermore, the neutral EWIMP should has a charged 
SU(2)l partner, x ± . When the EWIMP is heavier than 
EW scale, their masses are almost degenerate, and the 
unsuppressed transition between the two-body states of 
2x and x~X + may play an important role in the 2%° 
pair annihilation. 

In this letter we reevaluate the pair annihilation cross 
sections of the EWIMP 's, for the two cases that the 
DM is a component of two SU(2)i-doublet fermions or 
of an SU(2)i-triplet fermion. These correspond to the 
Higgsino-like and Wino-like DM's in the SUSY models, 
respectively. Most interesting fact we find is that the an- 
nihilation cross sections to the gauge boson pairs for the 
SU(2)i-doublet (triplet) DM suffer from a zero energy 
resonance around m ~ 6(2) TeV, whose binding energy is 
zero [5] under the potential. Therefore, the cross sections 
would be enhanced significantly compared with ones in 
the perturbative estimations for 1(0.5) TeV. Fur- 
thermore, it is found that the cross section for 2%° — * 27, 
which is usually suppressed by a one- loop factor, becomes 
comparable to the other tree-level processes, such as 
2x° -> W + W~, around the resonance. This means that 
the mixing between the two-body states of x~X + an d 2x° 
is maximal under the potential. Due to the explosive en- 
hancement of the cross sections, the SU(2)i-triplet DM is 
already partially constrained by the EGRET observation 
of the 7 ray from the galactic center, and the future 7 
ray searches may have sensitivity to the heavier EWIMP 
DM. 

First, we summarize properties of the EWIMP DM's. 
If the DM has a vector coupling to the Z boson, the 
current bound from the direct DM searches through the 
spin-independent interaction [6] is stringent. This means 
that the EWIMP DM should be a Majorana fermion or 
a real scalar if it is relatively light. Here we consider a 
former case for simplicity. 
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A simple example for the EWIMP DM's is a neutral 
component of an SU(2)i-triplct fcrmion (T) whose hy- 
percharge is zero. This corresponds to the Wino-like 
LSP in the SUSY models. It is accompanied with the 
a charged fermion, x^- While x° and are almost 
degenerate in mass in the SU(2)l symmetric limit, the 
EW symmetry breaking by the Higgs field, h, generates 
the mass splitting, 8m. If 8m comes from the radiative 
correction, 8m ~ \j2a 2 (mw — c^mz) ~ 0.18 MeV for 
m ^> m\y and mz- Here, myy and mz are the W and Z 
boson masses, respectively, and Cw(= cos9 w ) is for the 
Weinberg angle. Effective higher-dimensional operators, 
such as h 4 T 2 /A 3 , also generate 8m, however they are 
suppressed by the new particle mass scale A. The ther- 
mal relic density of the DM with mass around 1.7 TeV is 
consistent to the WMAP data. 

Another example for the EWIMP DM's is a neutral 
component of a pair of SU(2) ^-doublet fermions (D and 
D') with the hypercharges ±1/2. This corresponds to 
the Higgsino-like LSP in the SUSY models. The x° is 
accompanied with a neutral Majorana fermion, x'°, m 
addition to a charged Dirac fermion, x ± - They are again 
degenerated in mass in the SU(2)l symmetric limit. The 
mass difference is generated by the effective operators, 
such as h 2 D 2 /A, via the EW symmetry breaking. The 
thermal relic density of the DM explains the WMAP data 
when the mass is around 0.6 TeV. 

In the current Universe the DM is expected to be 
highly non-relativistic as mentioned before. In this case, 
the perturbative pair annihilation cross sections of the 
EWIMP DM may have bad behaviors if the DM mass is 
heavier than the weak scale. One of the example is the 
annihilation cross section to 27 at the leading order. The 
process is induced at one-loop level, and the cross section 
is 4(l/4)7ra 2 a|/m^ for the SU(2) L -triplet (doublet) DM 
in the SU(2)l symmetric limit. The cross section is not 
suppressed by 1/m 2 , and the perturbative unitarity is 
violated when m is heavy enough. 

The NR effective theory is useful to evaluate the cross 
sections in the NR limit. In Rcf. [7] we studied the NR 
effective theory for the EWIMP in a perturbative way 
and found that the trouble in the cross section to 27 is 
related to the threshold singularity. In order to evaluate 
the cross section quantitatively, we have to calculate the 
cross section non-perturbatively using the NR effective 
theory [12]. 

For evaluation of the annihilation cross sections for 
heavy EWIMP, we need to solve the EWIMP wave func- 
tion under the EW potential. In this paper, we show the 
formulae for evaluating the cross sections in the SU(2)l- 
triplct DM case. Those for the SU(2)i-doublet case will 
be shown in the further publications [8]. 

The NR effective Lagrangian for two-body states, 
</>jv(r)(~ l/2x°x°) and <fc(r)(~ x~X + ), is given as 

C = i* T (r) ^E+^pjl- V(r) + 2iTS 3 (r^j *(r) , 

(1) 



where 3?(r) = (<frc (r), 4>n(t)), r is the relative coordinate 
(r = |r|), and E is the internal energy of the two-body 
state. The EW potential V(r) is 
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In this equation we keep only 28m in (1,1) components 
in o rder to calculate the DM annihilation rate up to 
0{y/5m/m) [7]. T is the absorptive part of the two-point 
functions. Note that a factor of 1/2 (l/\/2) is multiplied 
for V22 and ^2 (V12 and T12) since 4>n is a two-body 
state of identical particles. Thus, T22 (Tn) is the tree- 
level annihilation cross section multiplied by the relative 
velocity v and 1/2(1). Since the SU(2) L -triplet DM is 
assumed to be a Majorana fermion, the 1 ,5- wave gauge 
contribution to T is relevant to the NR annihilation, and 
then, 



(3) 



The annihilation cross section of x~X + or 2x° to the 
EW gauge boson pair can be expressed using the two-by- 
two Green function, G(r, r'), which is given by 

( { E + ^r) 1 ~ v(r) + 2iTS3 ^) G ( r ' r ') 

= (5 3 (r-r')l. (4) 

Due to the optical theorem, the long-distance (wave func- 
tion) and the short-distance (annihilation) effects can be 
factorized [9]. The annihilation cross sections to VV 
(U, V = W, Z, 7) are written as 




{av)vv = Ci^2 ab T ab \ vv , xA a Al , 



(5) 



where i represents the initial state (i = and ± for 
2%° and x~X + pair annihilation, respectively) and A a — 



fd 3 r 



r e- lkr (£+V 2 /M)G M (r, 0) with k = ^JmE = mv/2. 



Here Cq = 2 and c± = 1, where Cq is a factor needed to 
compensate the symmetric factor for T and V. r a (,|yy 
is the contribution to T a h from the final states VV'. 
It is clear that if the long-distance effect is negligible, 

(av)vv = CiTu\vv- 

The S'-wave annihilation is dominant in the NR annihi- 
lation. Thus, the Green function is reduced to G(r, r') = 
g(r,r')/rr' . Similar to the case in one-flavor system, we 
find that g(r, r')/rr' is expressed by the independent so- 
lutions of the homogeneous part of the Eq. (4), g>(r)/r 
and g<(r)/r, as 



l(r, r') 



47T 



g>(r)gl(r')6(r - r') 



m 



+^S<(r)gl(r')0(r'-r) . (6) 

The solutions g> (r) and g< (r) are also two-by-two matri- 
ces since $(r) has two degrees of freedom. The boundary 
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FIG. 1: The x°-P a i r annihilation cross sections to 27 and W + W~ when 5m = 0.1 GeV (slid lines) and 1 GeV (dashed lines). 
X° is the SU(2)_L-triplet or doublet DM. Here, v/c = 10 -3 . The leading-order cross sections in the perturbation are also shown 
for 8m = (dotted lines). 



conditions at r — are g<(r)| r _,o = 0, g<(r)| r _>o = 1, 
and g>{r)\ r ^o = 1. In the following, we assume E < 
25m so that a pair annihilation of x° does not produce 
on-shell x~X + ■ As the result, 
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die ikr d 2 e ikr 



(7) 



In this case, the x°~P a i r annihilation cross sections are 
(av)vy = CiJ2ab^ab\vv'd a dl, as expected. It is enough 
to calculate d in order to evaluate the cross sections. 

In Fig. (1) we show the annihilation cross sections of 
the SU(2) L -triplet DM pair to 27 and W+W~ as func- 
tions of to. We evaluated the cross section numerically. 
Here, we take v/c = 10~ 3 , which is the typical averaged 

The qualitative behavior of the cross sections around 
the first resonance may be understood by approximating 
the EW potential by a well potential. Taking cw = 1 for 
simplicity, the EW potential is approximated as 

-bi^/2a2mw 




V(r) 



25m — b\a.2mw 
—bi^/2a 2 mw 



(8) 



for r < R(= {b2mw)~ 1 )- Here, b\ and 62 are numerical 
constants. By comparing the annihilation cross sections 
to 27 in this potential and in the perturbative calculation 
for small to, we find b\ = 8/9 and b 2 = 2/3. Under 
this potential, two-body states 2\° and x~X + have the 
attractive and repulsive states, whose potential energies 
are A± = 1/2(V U ± y/V^ + 4V 2 2 ) with Vy(i, j = 1, 2) 
elements in V. The attractive state is — sin 9<fic +cos 9(f) n 
with tan 2 9 = A_/A + . 

When 5m <C bia2mw/2(~ 1 GeV), 9 is not suppressed 
by 5m and x~X + an d ^X° are mixed under the potential. 
In this case, the cross section to 27 is given as 



(0-1O27 



47TCr 



1 



1 



(9) 



9m 2 \cos(k-R) cosh(fc + i?) 

where k\ = |A±|m. Here, we neglect the T term contri- 
bution to the wave function for simplicity and take E ~ 0. 



velocity of the DM in our galaxy, and 5m = 0.1 GeV and 
1 GeV. The perturbative cross sections are also plotted. 
Large 5m leads to unreliable numerical calculation for 
large m, and then some curves are terminated at some 
points. However, 5m should be suppressed around the 
regions. 

When to is around 100 GeV, the cross sections to 27 
and W + W~ are almost the same as the perturbative 
ones. The cross section to 27 is suppressed by a loop fac- 
tor there. However, when to ^ 0.5 TeV, the cross sections 
are significantly enhanced and have the resonance struc- 
ture. Especially, the cross section to 27 becomes com- 
parable to that to W + W~ around the resonance. This 
suggests that the 2%° state is strongly mixed with x + X~ ■ 
The cross section (9) is reduced to 47ra 2 a|/TO 2 ^ for 



a 2 m^~ mw- On the other hand, it is not suppressed by a 
one-loop factor for a 2 m^ mw and has a correct behavior 
as ~ 1/to 2 in a heavy to limit. When k-R = (2n— 1)tt/2 
(n = 1, 2, • ■ ■ ), the zero energy resonance, whose binding 
energy is zero, appears and the cross section is enhanced 
significantly. In Fig. (1), the n-th zero energy resonance 
appears at to = to*-™' ~ n 2 x to/ 1 ), while the well poten- 
tial predicts mS n > ~ (2n — l) 2 x to/ 1 ). We guess that the 
Yukawa potential might be approximated better by the 
Coulomb potential for the higher zero energy resonances. 

When the zero energy resonance exits, the cross sec- 
tions av are proportional to v~ 2 for t < 1. However, 
this is not a signature for breakdown of the unitarity. 
We find from study in the one-flavor system under the 
well potential V that when v <C mVT, av is saturated 
by the finite width T and the unitarity is not broken. 

We also show the annihilation cross sections for the 
SU(2) L -doublet DM in Fig. (1). The SU(2) L -doublet DM 
has the smaller gauge charges compared with the SU(2) j,- 
triplet DM. As the result, the cross section is smaller, and 
the first zero energy resonance appears at 5 TeV. 

The enhancement for the DM annihilation rates gives 
significant impacts on the indirect searches for the DM 
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FIG. 2: a) Line 7 flux from the galactic center in cases of the SU(2)L-triplet and doublet DM's. We take the averaged velocity 
of the DM v/c = 1CT 3 , J = 500 and AS7 = 10" 3 . The leading-order cross sections in the perturbation are also shown (dotted 
lines), b) Contour plot of continuum 7 ray flux from the galactic center in the unit of cm~ 2 sec _1 GeV _1 . Here, we take 
8m = 0.1 GeV, and other parameters are same as in a). The shaded regions correspond to S/B > 1. 



in the cosmic ray. In the following we discuss the 7 ray 
search from the DM annihilation in the galactic center 
and the future prospects. 

The line 7 from the pair annihilation to 27 or Z7 at 
the galactic center is a robust signal for the DM. Also, 
W~W + and 2Z° final states produce the continuum 7 
spectrum through ir° — > 27, and the observation may 
constrain the EWIMP DM. The 7 flux, * 7 (.E), is given 
as 



dE 



9.3 x lO-^cm^sec^GeV- 1 x JAil 



/lOOGeVV 



,dN 



vv' 



\ m 



J ^ dE 

7 vv 



'vv 



10 _27 cm 3 sec -1 



.(10) 



where N vv is the number of photons from the final state 
In Fig. (2a) we show the line 7 flux from the galactic 
center in cases of the SU(2)i-triplet and doublet DM's. 
Here, we take dm = 0.1, 1, 10 GeV. We also show the flux 
obtained by the leading-order calculation for comparison. 
The ACT detectors have high sensitivity for the TeV- 
scale 7 ray. MAGIC and VERITAS in the north hemi- 
sphere might reach to 10 _14 cm _2 s _1 at the TeV scale 



VV and (av) is the averaged cross section by the velocity 
distribution function. Note that the angular acceptance 
of the detector, Afl, is 10 -3 typically for the ACT detec- 
tors. The flux depends on the halo DM density profile p 
through 



JAfl 



1 



8.5kpc J 1-o.s 

An 



dfldl 



P 



,0.3 GeVcm- 3 



(11) 



where the integral is along the line of sight. J is stud- 
ied for various halo models, and 3^ J^IO 5 [4]. The 
cuspy structures in the halo density profiles, which are 
suggested by the iV-body simulations, tend to give larger 
numbers to J. In the following we take a moderate value 
as J = 500, which is typical for the NFW profile [10]. 
while CANGAROO III and HESS in the south hemi- 
sphere to 10 _13 cm _2 s~ 1 [4]. These ACT detectors are 
expected to cover the broad region. 



In Fig. (2b), the contour plot of the continuum 7 flux 
from the galactic center is presented. For dN vv /dE 
we use the fitting functions given in [4]. Shaded regions 
correspond to S/B > 1. In order to evaluate the back- 
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ground B, we assume a power low fall-off in the energy 
for the diffused 7 ray flux ^ BG {E) as d^ BG (E)/dE = 
9.1xlO- 5 cm- 2 sec- 1 GeV- 1 x(£;/l GcV)- 2 7 Arj [4]. The 
EGRET experiment has observed the diffused 7 ray emis- 
sion from the galactic center up to about 10 GeV [11]. 
Even the small regions around the resonances in addi- 
tion to the triplet DM with m x = 100 GeV are already 
constrained by the EGRET observation. The GLAST 
satellite detector, which will detect 7 ray with 1 GeV< 



E <300 GeV, have more sensitivity to the around the 
region. 
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